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I. INTRODUCTION 

Throughout this paper CT and DT refer to continuous time and discrete time respectively, units are taken as 
c = U = 1 and the metric tensor has components ?7 M „ = diag(l, — 1, — 1, — 1). The symbol — c denotes an equality 
holding modulo the equations of motion. In recent years there has been much interest in the discretization of time for 
a variety of reasons. These include the study of integrable systems Jl| , motivated by the question of the integrability 
of various dynamical models. Systems studied are often point particle models with special properties which make the 
process of discretization successful, such as the Toda lattice field theories [3-4], considered as approximations 
to CT field theories. Generally these tend to look like variants of lattice gauge theory, with the important difference 
that time is not Wick rotated; lattice gauge theory, which has as its objective the CT limit a — > 0, where a is 
the lattice spacing and which is regarded as the place where physics occurs. The value a ^ is taken only as a 
■ calculational device, suitable for providing a regularization procedure and for computer simulation to some given 
I/"") | order of numerical accuracy. In this scheme discretization is also applied to the spatial coordinates, and because time 
. is Euclidean, the theory uses a four dimensional Euclidean space. On a more fundamental level it has been speculated 
by many authors that perhaps space and/or time are not continuous on very short scales. In the history of quantum 
field theory, problems associated with operator products at the same spacetime point were sometimes addressed by 
"point-splitting" techniques, wherein two or more field operators at the same point were defined at separate points 
and only after the calculations was the separation taken to zero. This was a technique used for example by Schwinger 
i-pH ■ in his famous 2-dimensional model of QED ||. More recently there has been a strong suspicion that at the Planck 
scale the usual view of spacetime breaks down and novel new ideas should be contemplated. In a series of papers [3-6] 
attention was focussed on the consequences of discretization of time but not of space. Such an idea was discussed by 
't Hooft in 2 + 1 dimensions JlQ] , with the suggestion that under certain circumstances, there may be variants of the 
Regge calculus approach to General relativity |T1J wherein time becomes discrete but space does not. More recently, 
discrete time has been discussed as a direct consequence of quantum processes in action p2| . In this approach, time 
jumps whenever a factor state in the quantum state of the Universe undergoes a test resulting in information gain. 
These factor states can in principle include attributes such as momentum, so there is no necessity to discuss discrete 
space as well, although that is possible. In this article only temporal discretization will be discussed. In the system of 
DT mechanics discussed here, it is considered an exact form of mechanics, with its own consistent laws of motion and 
dynamical invariants, and not only in the CT limit. This is a fundamental difference between this DT mechanics and 
conventional lattice gauge theory. All of the operators in our DT quantum theory have to be good, in the language of 
lattice gauge theory. It is the case that any conventional continuous time theory can be discretized, simply by looking 
at the system at discrete times. For example, integrating the Lagrangian between chosen times (temporal nodes) gives 
Hamilton's principal function (the Hamilton- Jacobi function) , which can be regarded as a convenient discretization of 
time, because it depends only on what the system is doing (the co-ordinates) at the end points of the temporal interval. 
Our original motivation for discretizing time arose when computer simulation of a soliton theory was undertaken [ p"3"| . 
In CT theory, various integrals such as total mass, linear and angular momentum and charge remain invariants under 
dynamical evolution. However, it was soon found that a naive discretization of the equations of motion did not lead 
to the naive discretizations of these integrals remaining invariant during the simulations. This has nothing to do with 
numerical inaccuracy, but with the principles being employed. It was perceived that discretization should perhaps be 
done in a more carefully controlled, principled way, so that a genuine discrete dynamics with true invariants emerged. 
In this paper some of the methodology discussed in j?|] is applied to a well-known field theoretic model with an 
underlying non-abelian group structure, the Skyrme model. This is not a trivial application, for it is not obvious 
at all that every model in CT field theory carries over easily into a DT analogue theory with DT analogues of the 
CT structures remaining intact. An important difference is that there may be no obvious DT analogue of the CT 
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Hamiltonian. It is the merit of integrable systems that such an analogue should be found, but in general this will not 
be so. The obvious reason is that by definition, there is no continuous translation in time in DT mechanics, so there 
is no DT analogue of the CT generator of such a transformation. Fortunately, the basic idea of a Noether invariant 
survives in DT mechanics. If there is a continuous symmetry of the system function (the DT analogue of the CT 
Lagrangian) then there will be a corresponding invariant, referred to here as a Maeda-Noether invariant [[uj. In the 
DT Skyrme model, the symmetry which generates the axial and vector charges survives and the corresponding DT 
charges are discussed here. We shall also discuss T.D. Lee's version of DT mechanics, which automatically gives a 
D T analogue of energy. 



II. DT MECHANICS 



In the models discussed here the dynamical variables are defined only at instants of time nT where n is an integer 
and T is some fixed time interval. In these models there are no gauge fields defined on the temporal links connecting 
successive instants of time, so the formalism presented here takes on a simpler form than would be the case in say 
QED or QCD. The dynamical variables here are real valued fields ip a , a = 0,1,2,3 with suitable constraints. The 
index a is not a Lorentz index, but it is lowered with the metric tensor components rj^. The DT analogue of the CT 
Lagrange density C = C (ip a , d^ip 01 ) is called the system function density and takes the form 

S" = S«,^ +1 ,V<,V< +1 ). (1) 

It has the physical dimensions of an action density, not an energy density. The integral S n = J dxS n is called 
the system function, and although used as the DT analogue of the Lagrangian L in CT mechanics, is more like a 
Hamilton's principal function. Applying Hamilton's principle to the action sum 

JV-l 



a mn [r] = J2 sn ( 2 ) 



n=M 



gives the second order equations of motion M 

n£= c n£, M<n<N-l (3) 

where 

Given an infinitesimal transformation of the fields tp a — > ip a + 5tp a then 

SS n = c S^ +1 K +1 - 6^K + V • a, (5) 

where 

F) ci 9 n 

A symmetry of the system function is a transformation of the fields which leaves the system function unchanged. 
Hence for such a symmetry 

<5< +1 H£ +1 - S^K + V • Sa n = 0, 
=> £< +1 n£ +1 -<5<n£+V-&r„= c 0, (7) 

which can be written in the forwards form 

D+S Pn + V • 5j n = c (8) 
using the equations of motion (||) . Here the forwards densities are defined by 

5p n = 5<p%n%, 5- Sn =T- l 8a n , (9) 
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with D+ = (U n — 1) /T , where U n is the temporal step operator defined by U n f n = f n +i- From (||) the corresponding 
dynamical invariant is given by 



sc n ^ J 8<t%n$, (io) 



assuming that the fields fall off sufficiently rapidly at spatial infinity. An alternative form is given by noting that 

S^+iK+i ~ S<PnK + V • Sa n =o 0, (11) 
giving the backwards DT equation of continuity 

D-8p n + V • S] n = c (12) 
where D~ = (l — U^ 1 ) /T and the backwards densities are defined by 

Spn = 5<pZ +1 TL% +1 , 5i n =T- l 8a n . (13) 

III. THE SKYRME MODEL 

The Skyrme model fl5[ | has soliton solutions with a number of properties suggestive of baryon physics. Its basic 
dynamical degrees of freedom are space-time fields U (x) which take values in 5/7(2). Because this implies constraints, 
it is often convenient to parametrize these fields in terms of an unconstrained isotopic triplet of real scalar fields 

7T = (7T 1 , 7T 2 , 7T 3 ) I 

T • 7T 

U = exp {it-it} — cos (| ir\) + i sin (| 7r|) , (14) 

\ir\ 

where r = (t 1 ,t 2 ,t 3 ) are the Pauli matrices. With the definitions 

U„ = d^U, L M = U+d^U = -d^U+U (15) 
the Lagrange density may be written in the form 



C s = -^TrL^ - ^Tr [L^L V ] [L»,L V ] (16) 



where F n is the pion coupling constant and the second term is known as the Skyrme term. In terms of the U fields 
this is equivalent to 

Cs = lf Tr U t U " ~ lh Tr { U Z U » U+>lUV - U+U U U +V U^} . (17) 

In addition to the standard Poincare symmetries an important symmetry of this Lagrangian is invariance under 
separate left and right SU(2) transformations; 

U^U' = AUB+ (18) 

where A and B are spacetime independent elements of SU(2). This generates the so-called axial and vector charges. 
In Appendix A the quaternion approach to the parametrization of the U variables is given. With this, U may be 
written in the form 

U = q^ = <p° + iT i ip i , (19) 

where qo = I2, the 2x2 identity matrix and qi = it 1 , i = 1, 2, 3 have the properties of the quaternions i, j, k. The 
four real fields are read off from (|l4| ) to be 

tp° ee cos (| tt|), (p* = sin (| 7r|) n • n = 1. (20) 
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Since there are only three independent parameters describing the elements of SU (2), the four components <pP are 
constrained to the surface of S 3 , the unit sphere in four dimensions, i.e. 

(^V = 1. (21) 

With this reparametrisation the Lagrange density becomes 

£s = yC^ Q dV" " ^d^d^ {dVdV Q " 0V0V} + \^{v a V a ~ 1) (22) 

where a 2 = bF%, 1 = e~ 2 and the Lagrange multiplier /i enforces the S 3 constraint rt|jh . Then the conjugate 
momenta are given by 

bp 

** = ^ = (23) 

where 

M a p = (a 2 - /3 2 di^di^) S af} + /3 2 di<p a di<p . (24) 
The constraints turn out to be second class in the terminology of Dirac and given by 

XL = <f a <p a -1«0, X2 = ^ Q vr a «0 (25) 
Then the non-zero Dirac brackets are evaluated to be 

{<,4) DB = (Kd - i£<p2) S 3 (x - y) . (26) 
It is these which should be used in the quantization of the fields. 



IV. THE SU(2) PARTICLE 
A. continuous time 



In this section the most basic variant of the Skyrme model is considered, which is to drop the Skyrme term and the 
spatial degrees of freedom. Then the Lagrange density reduces to the Lagrangian 

L = \a 2 Tr U + U, (27) 

where U = U (t) is a time dependent element of SU(2) and a — ^F n . The number of independent real dynamical 
variables is three and there are two alternative formulations: 



1. the it fields 



The U fields may be parametrized using three unconstrained real fields: U (t) = exp {it-it (t)} , where n (t) = 
F (t) n (t) is an element of 1Z 3 and n (t) is a unit 3— vector. Then 

U{t) = cosF + isinFr-n. (28) 

The mapping from the 1Z 3 space of the parameters 7rto SU (2) is many to one, with the vectors (F + 2kir) n, k an 
integer, mapping into the same point of SU(2). The Lagrangian ( |27| ) then reduces to 

L = \a 2 \P 2 + n-nsin 2 f\ + \[i (n n - 1) , (29) 

where a Lagrange multiplier is included to enforce the normalization condition on the unit vector n. The equations 
of motion are 
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F - sin F cos F n • n= c 0, (30) 
a 2 sin 2 Fh + 2a 2 sinFcosFFn= c fin, nn= c l. (31) 

In phase space the system has two second class constraints in the language of Dirac |L6) . Now define p, p to be the 
momenta conjugate to F and n respectively. Then these constraints take the form 

Xi=n-n-l«l, x 2 = n-p«0. (32) 

Following Dirac |L6| , the Dirac brackets can be constructed in the standard way giving the non-zero brackets 

{p,F} D = -l, {p l y} D = -8 l , J +n l n\ {p\p>} D = pV - p> V. (33) 

The total Hamiltonian is given by 

2ar 2ar sin b 

which is an invariant of the motion. Two additional invariants of the motion can be found using Noether's theorem 
by observing that the transformation 

U -»• U' = AUB+ (35) 

is a symmetry of the Lagrangian, where A and B are space and time independent elements of SU (2). Writing 
A ~ 1 + it ■ a, B ~ 1 + it ■ h, where a and b are infinitesimal then 

£F = (a - b) • n, <5n = cotF {a - b - n- (a - b) n} + nx (a + b) (36) 

to lowest order in the infinitesimal parameters. Now an application of Noether's theorem gives the conserved left and 
right charges 

L = n x p— pn — cot Fp, R = n x p+pn + cot Fp (37) 
in phase-space. In configuration space they take the form 

L = a 2 jsin 2 Fn x h Fn - cos F sin Fn j = A — V 

R = a 2 | sin 2 Fn x n+Fn + cos F sin Fnj = A + V (38) 

where 

A = a 2 sin 2 Fnxn, V = a 2 ^Fn + cosFsinFn^ , 
are conserved separately. These are known conventionally as the vector and axial charges respectively. 



The Lagrangian takes the form 



and gives equations of motion 



i.e. 



2. The tp fields 



L^ia 2 ^V* + i/i(^V Q -l) (39) 



a 2 Cp a = cm a , <^V*=cl, (40) 



^= c (<pW Q . (41) 

The Lagrangian is invariant to the global SU(2) transformation U' = AUB + where A and B are time independent 
elements of SU (2) . Now suppose 
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A ~ 1 + it ■ a, B ~ 1 + ir ■ b (42) 

where a and b are infinitesimal, then 

5ip° = (b - a) • ip, 5<p = (a - b) ip° + (p x (a + b) (43) 
Hence the conserved charge is given by 

Sp = (a - b) • V+ (a + b) • A 

where 

V = a 2 (tp°tp - tp°ip) , A = a 2 ip x ip (44) 

with V = c A = c 0. 

B. temporal discretization 

Turning to the temporal discretization of the SU(2) particle system, the problem reduces to choosing a suitable 
virtual path between temporal notes |(| . It was found that appropriate paths were of the form 

U nX = XU n+1 + (1 - X)U n (45) 

where the parameter A € [0, 1] interpolates temporal nodes and 

U n =q a ^, « = 1. (46) 

Given the Lagrangian 

L=\a 2 TrU+U (47) 

then the system function becomes 

S n = \Ta 2 D+^ n D+^ n + i Mn T (<^>« - 1) + \^ n+l T - l) (48) 

The equations of motion are 

< +1 -2^ + <_! = cM „T 2 <, «=d (49) 

This is equivalent to 

< +1 + <_! =c ^ + < (50) 

There is invariance under 

[/„ = A(7„B+, A, B £ 5(7(2) (51) 

Then if 

A~l+ir-a B~l + ir-h (52) 

then 

<fy?° = - (a - b) • <p n , 5<p n = <pP n (a-b) + <p n X (a + b) (53) 

Then 

£p„ = 5<n« = (a - b) • V" + (a + b) ■ A", (54) 

giving the conserved charges 



G 



A" EE a 2 D+(f n X if n = —ipn+i x ip n 

2 

V" ee a 2 [<p° n D+<p n - D+<p° n tp n ] = ?L (yl<p n+1 - ip° n+1 <p n ) (55) 

If these charges are non-zero, then there is a natural time independent frame in isospace given by the directions 
(V™, A", A™ x V"). The following argument simplifies the equations of motion and establishes the existence of an 
infinite hierarchy of quadratic invariants. First note that, regardless of the equations of motion, the quantity U^ +1 U n 
is an element of SU (2) , and so may be written in the form 

U+ +1 U n = gx^ (56) 

using the quaternionic notation discussed in Appendix A, where 

< = c a0 X ^ + l^n (57) 

and 

« = !• (58) 

In detail, the components of 4> A turn out to be 

< = *„ = v"-a", (59) 

where 

y „ = Z_yn = _ ^0 a „ = ]_ A n = x (60) 

or or 
^From this it follows that (^°) 2 is an invariant of the motion, namely 

($°) 2 = c ($°_ 1 ) 2 , (61) 

so we may write 

<P>Z+i=cCe n , (62) 

for some real constant C and where e n = ±1. It turns out that 

C 2 +vl+a.l = l, (63) 

which means that — 1 < C < 1. Hence the equation of motion can be written in the form 

<+i + ^-i=cC (e„ + e„_i)y#. (64) 

where the e n are of magnitude +1 but otherwise arbitrary. This arbitrariness can be traced to the use of the Lagrange 
multipliers fj, n in the system function (^|) and is not a feature that exists in the CT limit T — > 0. In the special 
case that e n = +1 Vn then the equation is recognized to equivalent to the DT harmonic oscillator discussed in ||. 
Moreover, the bounds on the constant C mean that the motion is never hyperbolic. In this case it is found that 

VnVn = 1, ^n+l^n = ( Pn+2 l Pn = 2 ^ ~ 1 

and generalizing this result gives the infinite set of invariants 

Vn+rnVZ = T m (C) 



where T m is a Chebyshef polynomial of Type 1 |17| In terms of the F, n description, the parametrization is given by 

U n = c n + is n n n (65) 
where c„ = cos F n , s n = smF n , n n • n n = 1, and then 

a" = s n s n+1 n n x n n+ i, v" = s n c„ +:L n„ - s„ + ic„n n+ i. (66) 
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V. THE a MODEL 



A. continuous time 

The model is now extended to include spatial dependence, but not the quartic terms in the original Lagrangian. 
We shall call this the a model. The CT Lagrange density is now 

C = \a 2 Trd^U + d^U, (67) 

which is equivalent to 

C = \a 2 d^ a d^ a + \n (^> Q - 1) . (68) 

The equations of motion are 

a 2 Uip a = c [up 01 , <p a ip a = c 1 (69) 

which reduce to 

□^= C (^D/)^. (70) 
The conserved energy-momentum tensor density is 

T» v = a 2 id' l ip ol d v ip ot - iT/^c^dV} ■ (71) 

The invariance of the Lagrange density to the same transformation as before allows the vector and axial currents to 
be determined. Under the infinitesimal transformation 

U -► U' = (1 + ir-a)U (1 - fr-b) (72) 

then the fields change according to the rule 

Sip = (b - a) • ip, 5<p = (a - b) cp° + <p x (a + b) , (73) 

giving the conserved currents 

= tp°&*<p - 9VVj ai " = 9A V x V- (74) 



B. discrete time 

The system function density for the cr-model is taken to be 

S n = \Ta 2 {D+v a n D+^ n - §V^ • Vipl ~ Vytf+i • V^ +1 } 

+ il>„ (« - 1) + \T» n+1 - 1) , (75) 

which gives equation of motion 

^K^^^ + ^J-TaV^^M, «= c l. (76) 

There is invariance of the system function density under the infinitesimal transformation flSSj ) , and this gives two 
conserved currents: 

D+V*+ftV<= c 0, DtK + diK=cO, (77) 

where 
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l.C 



A° = n„ x tp n 
A!, = 4 <, — — x ^„ 



<Pn+l 



8S r ' 



ddi<p° 



n+1 



8S n 



+ fn+1 



dS r - 



ddilf n+1 



T 1 ddnpr 



8S r 



ddiip 



n+1 



(78) 



V° ee ^ {^i - Wn+i} + |Ta 2 {^V 9 ^ - ^VVn} 

op 1 



A l n = -\a 2 diip n x <f n - \a 2 diip n+ i X <£>>i+l- 



(79) 



VI. DISCRETIZATION OF THE FULL MODEL 



We may rewrite the full Skyrme Lagrange density equation (22) in the form 

C s = \^ a M^ (di^) <pP-W (flj^) + |/x (v Q < - 1) , (80) 

where M af3 (8^) is given by equation (24) and the potential function W (ditp^) contains quadratic and quartic terms 
in the derivatives of the fields. A serious problem may arise if any term in a Lagrangian is greater than quadratic in 
the dynamical variables. It is possible in such a case that certain discretizations leads to implicit equations of motion 
which cannot be solved directly to give the future values of the variables from a knowledge of their previous values. 
This remark applies to a classical theory. It is interesting that quantization may improve the situation. Suppose we 
had a system with dynamical variable x n , and the DT equation of motion was an implicit one of the form 

$ (x n -x,x n ,x n+1 ) = c 0, n = l,2, ... (81) 

where <!> is some function. If we could not solve this equation, viz, find a unique expression of the form 

x n+ i = c ft (x n -i,x n ) , (82) 

then the implication is that there exists more than one solution consistent with the initial conditions. Alternatively, 
given £ n _i and x n +x, there would not be a unique DT trajectory (i.e., unique value of x n ) connecting these points. 
This is of course unsettling and contrary to classical principles, but not so in quantum theory. We just have to recall 
Feynman's path integral, which explicitly requires us to consider all possible trajectories between initial and final times, 
including non-classical ones. It will be seen from this that DT mechanics should accommodate quantum mechanics 
better than classical mechanics, except in simple cases such as the DTharmonic oscillator. The DTFeynman rules for 
a ip 3, scalar field theory was discussed in B, and it was found that the momentum space rules had softened vertices, 
due to the effects of temporal discretization, which is a potentially useful result. Analogous effects are expected from 
non-commutative spacetime theories. There are several possible discretizations of the Lagrange density d8C| ) , such as 

S n = \TD+ V ^M^ (i^ + D+<f% - TW (Idutt + 

+ i7> - 1) + il> - 1) (83) 

or 

S n = \TD+ V a n {M** + M al3 (d^+i) } D+rf 

-±T {W +W (dup®} 

+\T» (« - 1) + iT/i - 1) . (84) 
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Whichever form is chosen, the (implicit) equation of motion will be given by formula (3) . In any case, the symmetries 
discussed in previous sections will hold, because the variations are global. The conserved DT vector and axial charges 
can be readily worked out, and left as an exercise. There is now a guarantee that these are dynamical invariants, even 
though it may not be possible to put the equation of motion into explicit form. 

VII. DISCRETE TIME ENERGY 

The methods discussed above permit the construction of those invariants associated with continuous symmetries, 
such as linear and angular momentum, and various charges, but not energy. Fortunately, there is a way of extending 
the formalism to generate a DT analogue of energy. This was discussed by T.D. Lee p8[ . The method is to take the 
time intervals T n = i n +i — t n to be dynamical, i.e., subject to their own dynamical equations of motion. The first 
thing is to change the system function density according to the rule 

S n ^S «, V<, V< +1 ,T) -> S n (T n ) ee S «, V<, V< +1 ,T„) , (85) 

where T n is the temporal measure between nodes n and n + 1. The action integral now becomes 

AT-l /N-l \ 

Amn [T] = ]T S n (T n ) - A ]T T„ - T MN , (86) 

n=M \n=M ) 

where the Lagrange multiplier enforces the constraint that the sum of the temporal measures adds up to the fixed, 
total time Tmn- The equations of motion for the ordinary dynamical variables follows exactly the same pattern as in 
the regular case (constant T), but now there are the additional equations 



dT n 

This gives a guarantee that the object 



d S n (T n )= c \, n = M,M + l,...,N-l. (87) 



C" = ^S n (T n ) (88) 



is a dynamical invariant. Actually, the correct way to see what is happening is to see the equality of the C™ as 
dynamical equations which the T n must satisfy. In general, numerical analysis would be required to solve these 
extended equations, and this may be non-trivial. Formally, however, the problem of energy drift is completely solved 
by this technique. Finally, this approach does not cause a problem with other dynamical invariants, such as the vector 
and axial charges. 

VIII. QUANTIZATION 

Quantization hi DT particle and field theories can be readily discussed using the Schwinger action principle, and 
it should be possible to extend this method to the Skyrme Lagrangian discussed here. This will be left for another 
occasion. 
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x. APPENDIX A : QUATERNIONS 

Given the Pauli algebra 
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where I2 is the 2x2 identity matrix, define the quaternionic symbols 

q l = ~ir\ q°=I 2 . (89) 
These satisfy the quaternionic multiplication rule 

= -S m °+e ijk q k . (90) 

Now define upper and lower quaternion indices with conjugation as follows: 

q ti =(q°,q i ), q^ = {q ,qi) 

qo = q° = (««)*. (91) 

Then conjugation is equivalent to the following quaternionic index raising and lowering action: 

(<zT = <^ M* = <f- (92) 

The product rule (^) may be written in the compact form 

<?<t = c»\q x (93) 

where 

C q — 1, G g C q U, C q — OlJ 

"", »• C 0l ,=c%=%, c^ = e iifc . (94) 



c 



Note that c M q = rf v = rj^, where rf v are the components of the Lorentz metric tensor. The product of three 
quaternions is given by 



tfff = c^c^ (95) 

and similarly for higher powers of the quaternions. The Pauli matrices are 2x2 matrices and this permits us to define 
a linear mapping Tr, called the trace, on the quaternions to the reals defined by 



Trq» = 25%, Trq^q v = 2rf v 
Trq a q l = 0, Trq i q i = -2%. (96) 
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The trace operation and the product rule (|93|) permits all SU(2) expressions to be simplified in terms of the 
quaternions rather than 2x2 matrices. 
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